Interactions in Chaotic Nanoparticles: 
Fluctuations in Coulomb Blockade Peak Spacings 
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We use random matrix models to investigate the ground state energy of electrons confined to a 
nanoparticle. Our expression for the energy includes the charging effect, the single-particle energies, 
and the residual screened interactions treated in Hartree-Fock. This model is applicable to chaotic 
quantum dots or nanoparticles — in these systems the single-particle statistics follows random matrix 
theory at energy scales less than the Thouless energy. We find the distribution of Coulomb blockade 
peak spacings first for a large dot in which the residual interactions can be taken constant: the 
spacing fluctuations are of order the mean level separation A. Corrections to this limit are studied 
using the small parameter l/kpL: both the residual interactions and the effect of the changing 
confinement on the single-particle levels produce fluctuations of order A/y^kpL. The distributions 
we find are significantly more like the experimental results than the simple constant interaction 
model. 

PACS numbers: 73.23.Hk, 05.45. Mt, 71.10.Ca, 73.20.Dx 
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I. INTRODUCTION 



The ground state properties of electrons in nanoparti- 
cles or quantum dots — finite systems of fermions inter- 
acting via Coulomb forces — remain incompletely under- 
stood. One valuable way to probe these systems is via 
electrical transport (for reviews see Refs. |l|-|). The 
dominant effect is the suppression of the conductance G 
because of the charging energy needed to add an electron 
to the nanoparticle and so allow current to flow. This 
Coulomb blockade (CB) is an essentially classical effect, 
and the charging energy is simply Ec — e^/2C where C 
is the capacitance of the nanoparticle. The blockade can 
be lifted by adjusting the potential on the particle — with 
a gate voltage Vg , for instance — so that the energy for N 
electrons is the same as that for iV-|-l. The conductance 
through the nanoparticle has peaks at these degeneracy 
points; in fact, G{Vn) is periodic with a spacing of e^/C 
between the peaks j^,^ . 

At low temperature the electron remains coherent in 
the dot and so quantum interference effects must be in- 
cluded in the Coulomb blockade 1^-^. Two main effects 
occur: (1) the peak heights fluctuate from peak to peak 
because of varying coupling between the wavefunction 
in the nanoparticle and the leads, and (2) the spacings 
between the peaks fluctuate because of interference con- 
tributions to the ground state energy. Peak heights have 
been studied first theoretically [p|-p^ and then experi- 
mentally |19|-pl[ ; agreement is good once effects of tem- 
perature 19| and classical dynamics [l^[l^] are taken 
into account. 

In this paper we focus on quantum effects in the peak 
spacings. These are directly related to the interference 
contributions to the ground state energy. One obvious 
contribution is the single-particle energy levels caused by 



the confinement in the quantum dot. But there may also 
be contributions coming from the interactions in the dot. 
We will show that for typical parameters characterizing 
experimental quantum dots or nanoparticles the contri- 
butions coming from the interactions must be included. 

The simplest treatment of Coulomb blockade peak 
spacing fluctuations results from writing the ground state 
energy as the sum of the classical electrostatic energy 
and the energies of the occupied single particle states 
§ . Such an expression for the ground state energy of a 
many Fermion system has a long and venerable history 
in atomic and nuclear physics |2^-gJ]. Using £ for the 
energy in the presence of the gate voltage while reserving 
E for = 0, we have for N particles 



£N{Vg) ^En- {Cg/C)eNVg 
En = {e^/2C)N^ + ^ 



(1) 
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where Cg is the capacitance of the gate to the dot and 
{eio-} are the single-particle energy levels. This treat- 
ment, known as the constant-interaction model (CI), is 
still an important point of reference in the field. The gate 
voltage V* at which a peak in the conductance occurs is 
given by the condition £M{Vg)^£N+i{Vg). The spacing 
between two peaks is, then, proportional to the (discrete) 
inverse compressibility A 

{V;)n^n+i - {V;)n^i^n (X (2) 
En+1 + En^i — 2En = An ■ 

In the case of the CI model, one gets the simple predic- 
tion I25II 



An = eVC 

An = e'^/C + {eN/2+1 - eAr/2) 



for odd TV 

for even N . (3) 
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Note the drastic odd/even effect in the quantum correc- 
tion to the classical result e^/C because of the spin of 
the electron. 

A number of experiments have studied the distribu- 
tion of CB peak spacings in semiconductor quantum dots 
p5|-p9[, a particularly flexible and convenient type of 
nanoparticle. The overall result has been striking dis- 
agreement with the CI model. The main discrepancy is 
that there is no sign of the strong odd/even effect pre- 
dicted by Eq. ( p[), though there is evidence of a weak 
odd/even effect ]29|p0| . In addition, in some of the ex- 
periments the widths of the observed distributions are 
significantly larger than what is predicted using the CI 
model and assuming the single-particle energy levels 
are distributed according to the appropriate Gaussian 
random matrix ensemble. These discrepancies indicate 
unambiguously that interaction effects play a significant 
role in the fluctuating part of the Coulomb Blockade peak 
position. 

Theoretically accounting for this disagreement is difh- 
cult because of two intrinsic characteristics of the prob- 
lem: (1) the fairly large number of electrons contained in 
the dot (typically A^- 100-1000), and (2) the long range 
character of the Coulomb interaction which requires that 
screening has to be correctly taken into account. The first 
of these points makes it delicate to extend to realistic dots 
the results obtained by exact multiparticle diagonaliza- 
tion ||2^,^, which can only treat problems with a few 
electrons ( <,15). To reach the experimental number of 
particles in a numerical approach, one has to rely on some 
self-consistent technique. One approach is Hartree-Fock 
which has been recently used in quantum dots [p2|-p5[. 
Then, however, the second point conies into play. In- 
deed, screening is a higher-order interaction effect, and 
it is well known for the Coulomb interaction that the ex- 
change terms can be significantly over evaluated in self- 
consistent Hartree-Fock (for the bulk see Ref. but 
this applies also for restricted geometries). This makes 
it rather difficult to interpret numerical self-consistent 
Hartree-Fock calculations using the bare Coulomb in- 
teraction as input. Another possible self-consistent ap- 
proach is, of course, density functional theory in the local 
density or local spin density approximation ||37|-|39|| . Here 
the screening is properly handled but approximations are 
made in the exchange term. Since, as we explain in de- 
tail in this paper, the subtle interplay between exchange, 
residual direct interactions, and single-particle level fluc- 
tuations is critical in describing quantum dot physics, 
results from this method are also difficult to generalize. 

An alternative theoretical approach is to formally carry 
out an RPA screening calculation and then use a short- 
range screened interaction in the subsequent evaluation 
of the magnitude of various effects |40|-|4^j8|. The na- 
ture of the screened interaction used is, then, important. 
Since experimentally the system size is much larger than 
the screening length, the most natural choice for the in- 



teraction is the bulk RPA screened interaction. For mo- 
mentum smaller than 2kp this is correctly approximated 
by the Thomas-Fermi expression. Throughout this paper 
we consider two dimensional systems, in which case 

^^TF(r) = j J^j2 ^TF(q) exp[iq-r] , 



Vtf (q) 



l + rs'{\q\/k,)/V2 



(4) 



where is the mean density of states (including spin de- 
generacy), and =ro/ao with Trrp the average area per 
electron and ag the Bohr radius in the semiconductor. 
It should be born in mind that although Vtf is propor- 
tional to Ts for small r^, it becomes largely independent 
of when >, 1, which is the case experimentally. In 
fact, within any RPA-like approach the interaction can- 
not become stronger than 



zero-range 



(r) = ,.-'6{r) 



(5) 



which we refer to as the zero-range limit of the potential. 
Once screening is included, the problem is that the inter- 
action is small, in a sense that we shall make more ex- 
plicit below, even in the strongest zero-range limit. Thus 
the self-consistent approach should give roughly the same 
answer as the CI model with first-order perturbation in 
the screened interaction added. This is exactly the pic- 
ture one has of weakly interacting Landau quasi-particles, 
with the obvious drawback that it is not straightforward 
to extract large contributions using this approach (see 
e.g. Refs. H^). 

In order to obtain a large interaction effect, one pos- 
sibility is to give up the Fermi-liquid-like approach and 
to assume, for instance, that for some reason screening 
is not as efficient in a nanoparticle as in the bulk, and so 
conclude that large zero-range interactions, not amenable 
to a perturbative treatment, should be used pSpl] , ^ , ^ . 
Before taking this drastic step, it seems reasonable, how- 
ever, to see how far one can go within the more standard 
Fermi liquid description. Indeed, even in this framework, 
we know that the CI model is inadequate. 

An accurate treatment of the spacing or inverse com- 
pressibility certainly requires that the average effect of 
the residual interactions be added to the CI model. By 
"residual interactions" we mean interactions beyond the 
simple classical effect taken into account in the con- 
stant interaction model and, in particular, interactions 
between the quasiparticles after the bare electrons are 
dressed through screening. While they make a smaller 
contribution to the total energy than those considered in 
Eq. (0), residual interactions make a contribution of or- 
der the mean level separation to the second difference A, 
comparable to the single-particle contribution in Eq. (|^) 
p0| , ^2| -|4^. In the limit of a large nanoparticle — one 
whose typical dimension L is many times the electron 
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wavelength, kpL ^ 1 — only the average effect of residual 
interactions needs to be added to the CI model p2|-^5|. 
We shall below refer to this approximation as the con- 
stant exchange model. The intuitive argument for this is 
that the effect of residual interactions involves integrat- 
ing wavefunctions over the entire volume of the nanopar- 
ticle and so self-averages: the particular characteristics of 
any given wavefunction are simply integrated out. Even 
if one neglects the fluctuations of the residual interac- 
tion terms, because of fluctuations of the single-particle 
levels, the "constant exchange" contribution can modify 
the total spin of the nanoparticle by favoring non-trivial 
occupation of the single-particle orbitals and, as a conse- 
quence, will significantly modify the fluctuations of the 
peak spacings The ground state spin of a large 

nanoparticle in this regime has been examined in detail 
1^,11 -47 1, but curiously the corresponding distribution 
of Coulomb blockade peak spacings has not appeared ex- 
plicitly (see however Ref. [^). Therefore, before turning 
to effects caused by fluctuations in the interactions — our 
main interest in this paper — we present below the distri- 
bution in the constant exchange limit (kpL 3> 1) to use 
as a reference point. 

In the experiments, however, kpL is not very large be- 
cause of the need to have the mean level separation larger 
than the temperature. This suggests that the variations 
in interactions caused by properties of individual wave- 
functions are important. The main goal of this paper 
is, therefore, to study the effect of fluctuations in the 
electron-electron interactions on the distribution of CB 
peak spacings. 

We limit our study here to the case of chaotic dynam- 
ics within the dot, for which a random matrix description 
of the non-interacting limit can be used |48j . Much as in 
the case of diffusive transport, completely chaotic dy- 
namics introduces an important additional energy scale 
(and only one): the Thouless energy, Etii = h/T, is the 
inverse of the transit time across the dot. It is gener- 
ally larger than the mean level separation A but smaller 
than the Fermi energy; we will assume A ^ E^h ^ Ep. 
On energy scales less than i?Th, a universal random ma- 
trix description holds while for larger energies (shorter 
times) the individual dynamics of the system comes into 
play. We will see that taking account of this additional 
energy scale is critical in describing CB peak spacings in 
nanoparticles. 

We focus on two effects in particular: (1) the fluctu- 
ation in the residual Hartree and Fock contributions to 
the ground state energy and (2) the change in the single- 
particle energies because of changes in the mean field 
potential as electrons are added to the dot [Q, an effect 
we refer to as "scrambling" . For two dimensional quan- 
tum dots, the changes in energy are evaluated to leading 
order in the smallness of the dot, the constant exchange 
model being used as the starting and reference point. 
We find that both effects contribute terms which are of 



order A/^/kpL, up to logarithmic corrections. The re- 
sulting distribution of CB peak spacings looks more like 
the experiment — it has a more Gaussian and symmetric 
form, for instance — but still deviates from the measure- 
ments in significant ways. 

The paper is organized as follows. Our starting point 
is a semiclassical expansion of the ground state energy 
described in Section O — an expansion in 1/kpL. The 



important large dot limit is presented in Section III 



The next two sections discuss the two main issues — 
fluctuations of the residual interaction contributions 
(Section |rv| ) and scrambhng (Section 0). Section VI 
takes a first step beyond the Gaussian model introduced 



in Section IV. Finally, our conclusions, including a dis- 
cussion of the experimental results, appear in Section 



VII 



II. APPROACH: SEMICLASSICAL 
CORRECTIONS 

In the Coulomb blockade through quantum dots we 
are faced with a classical theory which works remarkably 
well to which we want to add quantum effects to leading 
order. The small parameter is the standard semiclassical 
one, h. In our context, this is equivalent to an expansion 
in 1/kpL where kp is the Fermi wavevector and L is the 
typical size of the nanoparticle. Another useful way to 
view the corrections is in terms of the dimensionless con- 
ductance, g, given as the ratio of the transit rate hvp/L 
to the mean level separation A; in this case the expan- 
sion is in terms of 1/g, as in the diffusive mesoscopic 
regime For 2-dimensional systems, such as the ones 
we consider here, these two parameters are proportional 
(g = kpL/2Tr). 

We proceed by using the method of V. M. Struti- 
nsky [^,^ in which the dependence of many-body 
ground state quantities on particle number can be de- 
composed into an average and a fluctuating part. While 
the average part varies smoothly with particle number, 
the fluctuating part reflects the shell structure of the sys- 
tem. The smooth part is the bulk energy per unit vol- 
ume integrated over the finite-size system, and the oscil- 
lating contributions come from quantum interference ef- 
fects explicitly caused by the confinement. By supposing 
that the smooth part is known while the unknown oscil- 
latory contribution is a correction, Strutinsky introduced 
a physically motivated systematic approach to obtain the 
shell corrections |5^,|5^. This shell correction method is 
essentially a semiclassical approximation. It rests on the 
fact that the number of particles in the system considered 
is large, rather than on the interaction between the par- 
ticles being weak. (One must, of course, work in a regime 
where the smooth starting point is basically valid.) Since 
the quantum dots in which we are interested contain on 
the order of 100 electrons, they are a perfect place to 
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apply the Strutinsky method. We use the formulation 
recently developed specifically for quantum dots in Ref. 

H- 

Density functional theory guarantees that the ground 
state energy of a nanoparticle can be written as a func- 
tional of its density, E'dft = •^DFTi'^DFT]- Neglect- 
ing quantum interference effects in this functional cor- 
responds to a generalized Thomas-Fermi approximation 
(generalized because bulk local exchange and correlation 
can be included) which can also be written as a func- 
tional of a (smooth) density, i?GTF = -^gtf['t-gtf]- The 
quantum interference part of the energy, i^DFT — ^-otf, 
can be obtained approximately by expanding tidft about 
ncTF and solving the resulting density functional equa- 
tions order by order in the oscillating part of the density 
n-DFT — KGTF- To carry out this procedure an explicit 
form for the exchange-correlation functional is required. 
Based on the local density approximation, the result de- 
rived in Ref. ||5^ is 

£;dft ^ Egtf + AS^^^ + A£^^^^ (6) 
where the first and second order correction terms are 

N 

Ai?(i) = [Fcff [noTF]] = J2 

1=1 

A£;(2) = ly" drdr'n°'^r)Vsc{r,r')n°''%r') . (8) 

The first-order correction is the oscillatory part of the 
single-particle energy, fjp'^, calculated in the smooth 
Thomas- Fermi potential, Voff[nGTF]; that is, it is the 
sum of the deviations 6ei of the single-particle levels 
from their mean values. In the second-order correction, 
^osc ^YiQ deviation of the quantum mechanical den- 
sity calculated in the Thomas-Fermi potential from the 
smooth starting point, n°^'^ = n[T4ff [ncTF]] (r) — ?^gtf, 
and T4c(r,r') is the screened interaction in the nanopar- 
ticle. To the order that we are working, using h°^'^ here 
is equivalent to using tt-dft ~ '^gtf [ p2[ . 

In this approach, the ground state energy is, then, the 
sum of a classical contribution — the generalized-Thomas- 
Fermi result -Egtf — and two quantum contributions — 
A£;(i) and AE^^\ With the inclusion of only the first- 
order contribution, Eq. (^) is simply the CI model for 
the ground state energy. Because we have the correc- 
tions to this model, we can see exactly when it is appli- 
cable. The second-order correction, Eq. (^), has a nat- 
ural interpretation: the ripples in the density caused by 
interference (the Friedel oscillations in this setting) in- 
teract with each other via an interaction screened by the 
smooth Thomas-Fermi fluid already present. In fact, the 
bulk Thomas-Fermi interaction Eq. (^ can often be used 
for Vsc ■ Taking the — > cx) limit formally yields the zero- 
range interaction limit (^). This is equivalent to using a 
Hubbard Hamiltonian on a square lattice with hopping 
term t and interaction U with U/t = it. 



The result Eq. (^) is, however, inadequate for our pur- 
poses here because of the local density approximation 
used. In particular, the spin degree of freedom is han- 
dled poorly in this approximation: while Eq. ^ gives the 
"direct" interaction between the density ripples, an anal- 
ogous "exchange" term is missed. A better treatment 
of the interactions in density functional theory could 
presumably yield this additional term, by including off- 
diagonal terms in the density matrix, for instance. In 
fact, an RPA treatment of the interaction in the diffusive 
case yields an energy with exactly this structure p^,p[. 
Motivated by these physical considerations, we simply 
add the exchange term by hand to the expression for the 
energy i^DFT- 

For the purposes of deriving CB peak spacing distri- 
butions, we thus consider 3 contributions to the ground 
state energy: (1) the Thomas-Fermi energy, (2) the de- 
viation of the single-particle energy from its mean, and 
(3) a Hartree-Fock like treatment of the residual screened 
interaction. Without the residual interactions, a TV- 
particle state is, of course, a Slater determinant in the ba- 
sis of single-particle eigenstates ipi{r); it is characterized 
by the occupation numbers (noa, ni^, ■ ■ ■) where Ui^ = 
orl,(T = +lor— 1 labels the spin degree of freedom, and 

riicr = N . In this occupation number representation, 
the expression for the energy is, then, 

^tot({f^icr}) = ^GTF + ^ n^aSe., (9) 

i(7 

where 

M,, = J dvdv' \'4,,{vf Vtf{v - r') |V,(r')l' (10) 
N,, ^ J dvdv'Uv)r,{v)VMv - r')V,(r')V'*(r') (11) 

are the direct and exchange contributions, respectively. 
This is the starting point for our study of CB peak spac- 
ings. 

In order to derive the statistical properties of the 
CB peak spacings from Eq. (^), the statistical prop- 
erties of the single-particle eigenvalues and eigenfunc- 
tions must be known. For this purpose we assume that 
the single-particle classical dynamics in the nanoparti- 
cle is completely chaotic. In this case it is well-known 
that the single-particle quantum mechanics is accurately 
described by random matrix theory (RMT) on an en- 
ergy scale smaller than the Thouless energy i?Th (which, 
again, is the inverse time of flight across the system) 
p8| , p3[ . We consider only the case of no symmetry here, 
so that the energy levels {ei} obey Gaussian Unitary En- 
semble (CUE) statistics. Because of the spatial integrals 
in the expressions for and Nij , the correlations of the 
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wavefunctions are also needed. For the classical Gaussian 
ensembles, there is no spatial correlation for a given wave- 
function, except for that arising from the normalization. 
This is a direct consequence of the fact that no basis 
should play a particular role for these matrix ensembles. 
Here, however, because of the particular role played by 
the kinetic energy, the plane wave basis is special. More 
precisely, we should implement that each eigenfunction 
is localized in this basis on the energy scale i^xh- We 
shall therefore build our random matrix ensemble with 
the following requirements (a practical implementation 



is given in section VI) 



i) On energy scales smaller than the Thouless energy 
i?Th, one should recover standard GUE properties. 

ii) On scales larger than i?Th, an eigenstate of en- 
ergy e should appear localized on the energy surface 
|k| = y2rrie/fi. More precisely, we require that in 
a plane wave basis an eigenstate has an envelope of 
width 



5k = l/L 



where L is the typical size of the system. 



(12) 



iii) Finally, we have to implement that the fluctuations 
of energy levels, and also importantly of local den- 
sity, saturate at the Thouless energy. That is, if 
one considers a quantity A expressed as a sum over 
the energy levels in an energy window 5E about 
the Fermi energy, as the window widens the fluctu- 
ations of A increase until 5E sa i?Th beyond which 
point the fluctuations of A do not further increase. 

Except for the fact that the Fermi surface is given an 
explicit width « i?Th due to the finite size L of the sys- 
tem, items i) and ii) are very similar to Berry's modeling 
of chaotic wave functions which assumes that a ran- 
dom superposition of plane waves describes the correla- 
tion on scales smaller than the size of the system L. Item 
iii) is motivated by the known saturation of, for instance, 
fluctuations in the density of states of chaotic systems at 
an energy scale of order fi divided by the period of the 
shortest periodic orbit . Note that these requirements 
will yield a random matrix approach noticeably different 
from the ones discussed in Refs. |40 4^,^. 

In this model, the wave function statistics that we find 
are as follows: for a two dimensional dot of area A (the 
case that we consider throughout this paper) 



A(V'.(r)^*(r')) =<5,, Jo(fc,|r~r'|) 



(13) 



for |r — r'l <C L. (In a maximum entropy approach semi- 
classically restricted by the classical dynamics, this cor- 
responds to keeping only the direct trajectory between r 
and r'; the approach can be generalized to keep more 
classical paths [p5p^.) In addition, it can be shown 



that within our random matrix modeling, one has, up 
to l/{kpL) corrections 



(^,(ri)V;*(r2)V'.(r3)V*(r4)) = 
(V.(ri)V*(r2))(V'.(r3)^*(r4)) 
+ (^.(ri)V'*(r4))(V,(r3)V*(r2)) 



(14) 



Alternatively, these statistics can be obtained directly by 
assuming that the distribution of wavefunctions is Gaus- 
sian |55), though in this case special care must be taken 
with regard to the normalization constraint. 

There are alternate routes to our very natural starting 
point Eq. (H). One that has been developed recently is 
the "universal Hamiltonian" approach p6| , p7| , ^ J^ which 
uses RPA to treat the interactions, RMT properties of 
the single-particle wavefunctions, and the small param- 
eter 1/g to arrive at an effective Hamiltonian. Treating 
this effective Hamiltonian in the Hartree-Fock approxi- 
mation leads to the same expression for the ground state 
energy that we give above. 

We end this section by giving explicit expressions for 
the inverse compressibility, which is proportional to the 
peak spacing [Eq. (^], in two simple cases in which the 
standard up-down filling of the states is assumed. First 
consider the case of N even with n — N/2. For the stan- 
dard filling, it is a singlet state and the two neighboring 
states for iV— 1 and A^-l-l electrons are doublets. The 
peak spacing is then 

A'^ En even = A^Sgtf + Sen+i - Sen (15) 
1=1 

On the other hand if N is odd, still assuming the simplest 
up-down filling, the peak spacing is 



+ l,ri+l 



(16) 



The second difference of i?GTF is almost equal to the 
charging energy e^/C but there is a small correction from 
the residual interactions; this is discussed in detail below. 



III. LARGE DOT LIMIT 

The main simplifying feature in the large dot limit, 
kpL ^ 1, is that the fluctuations of the interactions can 
be neglected, much as for diffusive nanoparticles |5^-|59|] 
in the large g limit. The distribution of the spacings 
is determined, then, by the fluctuations of the single- 
particle levels and their interplay with the mean residual 
interactions. Of course, as has been studied previously 
p0| , p3| -|4^ , the ground state of the nanoparticle does not 
necessarily follow the simple filling used in Eqs. (p^- 
(p6|), and this effect must be included in finding the dis- 
tribution. 
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Let us consider the mean and variance of the My and 
Nij for levels near the Fermi level. We have, for instance, 
using Eq. ( p^ ) 

(M,,)^ y'drdr'(|V',(r)|^)yTF(r-r')(|^,(r')|') (17) 

+d,, I drdv'\{i;,{r)i;*{r'))f Vrpir - r') 

We will make the simplifying assumption that the 
nanoparticle is a billiard so that the average density is 
constant: (|'!/'(r)|^) — 1/A where A is the area of the par- 
ticle and the correlation function is given by Eq. ( [T^ ) . In 
this case the mean values are 



TF/fs 



(N.^,) - (18) 



where 

r2iT 



(T>TF)fs= / ^l>TF(fci.V2(l + COS0)) 



(19) 



var(M,^,) = J dridr2dr3dr4([|i^,(ri)|2 - A-^]VTF{ri - r^m^ir^)]^ - A'^] (21) 

x[|^,(r3)P - A-^]VTF{r3 - - A'^]) 

dridr2dr3duVTF{ri - r2)T/TF(r3 - r4) |(V'^(ri)^*(r3))n(V'j(r2)V':(r4))|' ■ (22) 



The first equality is due to the normalization of the wave- 
function, and the second derives from Eq. (p5|). Because 
of the decay of the autocorrelation of tp, Eq. (p^), the 
spatial integrals in ( ^l| ) are restricted to a region where 
Ti is close to r3 and likewise r2 to r4. These restrictions 
produce factors of the small parameter 1/kpL compared 
to the mean values. We will evaluate these factors be- 
low. However, for the large dot limit here, we see that to 
zeroth order in 1/kpL all the fluctuations in the residual 
interactions may be neglected |^^-^, as in the diffusive 
case [^6|-^ . Qualitatively, this is natural because of the 
averaging implicit in the integral in Mij and Nij: when 
the size of the system is much larger than the scale of os- 
cillation of '!/'(r), particular features of i/'i and ipj become 
less important compared to the mean behavior. 

For the CB peak spacing distribution, we thus con- 
sider the expression for the ground state energy (|) with 
constant (My) and {Nij). The problem is still compli- 
cated because of the interplay of the fluctuations in the 
single particle levels with these constants, which may, 
for instance, lead to = 1 or higher ground state spins 
pSj-^Sf. We thus evaluate the distribution of Coulomb 
blockade peak spacings numerically: Hamiltonians are 
taken at random in the GUE ensemble, and the lowest 
energy states for A^— 1, N, and iV+1 particles are ob- 



is the average on the Fermi surface of VTF(k — k'). In the 
zero-range interaction limit, these expressions simplify to 

{M,,)^{l + 6.,)A/2, {N,^,)^A/2. (20) 

The magnitude of the residual interaction effects is thus 
of order the mean level spacing A. (Note that the A 
that we use here and throughout this paper is the spac- 
ing of the orbital levels alone and so does not take into 
account the spin degeneracy factor.) Since energies of 
this order are critical in determining the spacing of CB 
peaks — note, for instance, the case of standard up-down 
filling Eq. ( p^ ) — the residual interaction terms must be 
included in any theory. In particular, the reason for the 
failure of the simple constant interaction model, Eqs. (|^) 
and (^), is now clear: it does not consistently keep all 
terms of order A. 

The variance or covariance of any of the My and Nij 
will involve the correlation function of a wavefunction. 
For example, consider the variance of Mi^j: 



tained by determining in each case the occupation num- 
bers {nio-} which minimize the expression Eq. (|^). 

The resulting distribution of CB peak spacings is 
shown in Fig. 1. We show it for two different strengths 
of interactions, the case of = 1 and the zero-range in- 
teraction limit. The first thing to note is that for both 
of these values the distribution differs substantially from 
that in the CI model, Eq. (1), with much less odd/even 
alternation present. Still, there is a (5-function in the 
distribution, as is clear in the cumulative distribution 
functions in panel (c), because the constant residual in- 
teraction term gives a rigid shift to the spacing for odd 
N and standard filling, Eq. (p6|). The origin of spacing 
in these plots is A^i^cTF + (M^^j). This is the mean 
spacing for adding a particle into the next available state 
neglecting all spin/filling effects; it thus corresponds to 
the classical charging energy e^/C. 

The fluctuations of the peak spacings are smaller in 
this constant exchange model than in the CI model. 
This is natural since the even and odd parts of the dis- 
tribution are brought closer together. The values are 
rms(A2£;Ar)/A = 0.58 in the CI model and 0.28 (0.24) in 
the present large dot case for = 1 (zero-range interac- 
tion) . 
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FIG. 1. The probability density of CB peak spacings in the large dot limit (fcpL — > oo). The total distribution (solid) as well 
as that for A'^ even (dashed) and odd (dotted) is given for two strengths of the interactions: (a) = 1 and (b) the zero-range 
interaction limit. The presence of a 5-function in the distribution is particularly clear in the cumulative distribution functions 
(the integral of the probability density) shown in (c). Insets show the probability of occurrence of ground state spins in the 
two cases. The spacing is in units of the mean level separation A, and the origin corresponds to the classical spacing f? jC . 



IV. FLUCTUATIONS OF RESIDUAL INTERACTIONS 

In order to evaluate the statistical effect of fluctuating residual interactions, we need to find the distribution of the 
Mij and Nij. We will show that in the limit of large but not infinite fc^L, the M's and iV's are Gaussian distributed 
and uncorrelated. 

The first step is to evaluate the variance of the M's and iV's. Using the random matrix ensemble introduced in 
Section ||, we find to a good approximation, 
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details for the case of var(Mi^j) are given in Appendix |^. The diagonal part of the direct residual interaction has an 
extra contribution because of the additional way in which the wavefunctions can be paired: 

(M ^ 9 rAf2 8 /-^^--/^ dq VTF(g)[VTF(0) + VTFiVm' ~ q')] 



In the zero-range interaction limit, the expressions for 
the variance of the M's and iV's simplify considerably, as 
for the mean values given in Eq. (^). In this case, we 
find 



var(My) = var(iVj^j) 



3A2 In(fcL) 
^ (fcL)2 



(1 + 3%). (27) 



Note that the decay of the wavefunction correlations ap- 
pearing in the variance produces a factor of 1/kpL in the 
rms compared to the mean. The In(fcL) factor is special 
to two dimensions; it comes from the 1/kL decay of the 
wavefunction correlator in this case. Thus in the large 
dot limit, kpL^l, the distributions of the M's and A^'s 
are narrowly peaked about their mean values. This is 
the justification for ignoring these fluctuations to a first 
approximation in the last section. 

Higher moments of the distribution involve coupling 



between more pairs of wavefunctions. The non-Gaussian 
part of the distribution is described by the cumulants, 
in which all of the wavefunctions are coupled in ways 
not present in the lower moments. In this case, the spa- 
tial integrals will have more restrictions coming from the 
wavefunction correlation ([l3|); these restrictions will in 
turn produce additional factors of the small parameter 
1/kpL. Thus in the large dot limit all higher moments 
can be neglected. In the case of third and fourth moments 
this can be easily verified by explicit calculation. Like- 
wise it is straightforward to check that cross correlations 
among the Af's and A^'s do not need to be considered, 
except for the obvious correlation between Mij and Nij 
which are equal in the zero-range interaction limit. We 
shall return to this point in Section VI. In this Section, 



however, we shall either take Mij ^Nij in the zero-range 
interaction limit, or neglect this correlation for rg — 1, 




FIG. 2. The probability density of CB peak spacings including the fluctuations in the residual interactions. The total 
distribution (solid) as well as that for A'^ even (dashed) and odd (dotted) is given for two sizes of quantum dot (N « 200, left 
column, or 80, right column) and two strengths of interaction (rs = l, top row, or the case of a zero-range interaction, bottom 
row). All of the features of the distribution are smoothed compared to Fig. 1. Note that the odd/even effect is quite clear in 
all four cases, and that the remnants of the 5-function are discernible. The spacing is in units of the mean level separation A, 
and the origin corresponds to the classical spacing e^/C. 
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and therefore take the My and iVy to be uncorrelated 
random Gaussian variables with mean and variance given 
above. 

We now use the statistics of the Mij, Nij, and sin- 
gle particle levels (assumed to follow the random matrix 
GUE ensemble as in Section III) to deduce the statistics 
of the CB peak spacings. It is clear from the expres- 
sions for the spacing in the simple up-down filling case, 
Eq. (p^, that sums such as ^i^l be important. 

Because of the saturation of the fluctuations, item iii) 
in the random matrix model of Section the variance 
of such a sum involves only levels within E^^ of each 
other. In calculating var(^^- Mij^, then, the sum should 
be taken over of order kpL independent M's: 



AY. 



M,, 



cx 



(28) 



The variance of other sums over the residual interaction 
terms will have the same dependence on kpL. Up to a 
•\/hi factor, we thus expect the contribution of the fluctua- 
tions of the residual interactions to the CB peak spacings 
to be of order A/^/kpL in rms. 

In order to evaluate the magnitude more accurately, we 
turn again to numerical evaluation. We consider a col- 
lection of i^Th/A levels with GUE energies and Gaussian 
residual interaction terms. We fill the levels with A^— 1, 
N, and N+1 particles successively, and find the occupa- 
tion numbers {ui^} minimizing the energy Eq. (^) in each 
case. Fig. 2 shows the resulting distribution of CB peak 
spacings. The fluctuations clearly act to smooth out all 
the sharp features of the fcj^L — > oo case: the (5- function 
becomes a finite peak, and the discontinuity at the origin 
is removed. The difference between the N even and odd 
cases is quite clear — even for the relatively small = 80 
(for which I/V^f-^— 0.21). The fluctuations of the resid- 
ual interaction terms lead, of course, to an increase in the 
fluctuations of the spacing; now rms(A^_Ejv)/A = 0.31 for 
rs — l and A^ = 80-200 compared with 0.28 in the absence 
of such fluctuations. While the inclusion of the fluctu- 
ations has clearly made the distributions more like the 
experiments, substantial differences remain. 



V. SCRAMBLING: ADDED ELECTRON 
CHANGES CONFINEMENT 

As electrons are added to the nanoparticle, the mean 
field potential which confines the electrons changes. This 
causes a change in the single particle levels, whose ener- 
gies appear in the expression for the ground state energy 
Eq. (||) , which then in turn change the CB peak spacing. 
Because we consider a chaotic quantum dot, the change 
in the energy levels will be unpredictable; hence, we call 
this effect "scrambling" . Note that it is an intrinsic effect 
connected to the charge of the electron not the geometry 



or environment of the dot. This effect has been studied 
previously for diffusive quantum dots in Ref. p2| ; here 
we give a derivation directly in the context of ballistic 
chaotic dots. 

We know from the Strutinsky approach (Section ||) 
that the change in the single particle levels should be 
found from the change in the smooth Thomas-Fermi po- 
tential, Eq. (0). So we first evaluate this change in poten- 
tial, and then use perturbation theory to find the change 
in the energy levels. Perturbation theory is justifled since 
the resulting shift is smaller than the mean level sepa- 
ration A by a factor depending on kpL. As a check, 
one can evaluate the magnitude of the fluctuation in the 
event of a complete scrambling of all the single particle 
levels. The properties of sums of random matrix eigen- 
values in Ref. |^0| imply a very large fluctuation of the 
peak spacings — of the order 2-4A which is larger than 
observed experimentally. 



A. Evaluation of S'VcS 

The generalized Thomas-Fermi problem is specified in 
terms of the density functional 

•^GTF [n] = Ttf [n] + £oxt ["-] + ^coui ["-] + £xc [n] (29) 
comprising four contributions: the kinetic energy 



TtfN = dr 



dn'e{n') 



(30) 



where e(n') is the maximum kinetic energy of non- 
interacting particles with density n' , the confinement po- 
tential defining the geometry of the dot 



(31) 



the Coulomb energy 

n(r)n(r') 



-'coul 



\n] — — / drdr' 
^ J 2 J 



(32) 



and the exchange-correlation potential £^c- The Thomas- 
Fermi density is found by solving 



6T, 



GTF 



6n 



.»^GTf] — fJ-TF 



(33) 



where n(r) is the electronic density and htf is found 
such that / n{r)dr — N. The effective potential is, then, 
defined by 

Vcs{r) = [ncTFl ■ (34) 

on 

Let SriQTF = '^■gtf(-^ + 1) — ?^gtf(-^) be the change in 
the Thomas-Fermi density when one electron is added to 
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the system. Bearing in mind that Eq. (|3^) is a classical- 
like equation for which SriQ^p is small, one can write 



dr 



[nGTF](r, r')(5nGTF(r') 



(35) 



rfr'^7^(r, r')(5r7,GTF(r') = S^tf 
on 



with SfjLTF such that J(5nGTF(r) = 1. In two dimen- 
sions, the Coulomb energy of the added charge will dom- 
inate any variations in the kinetic energy. One can 
therefore write the variation of the density as Jtt-gtf = 



^^GTF ~^ ^^GTF ~^ ■ 

trostatic problem 



with ^JT-GTF solution of the elec- 



dr' 



^1 I ^^/'^"'Gtf(i"') 



^(r,r')<5n^TF(r') 



r — r 



(36) 



B. Consequences for the peak spacing distribution 

The variation of the effective potential induces a 
change of the one particle energies. We let ejv denote 
the single particle energy of the electron, whatever 
orbital it is in (for N even and simple up/down fill- 
ing, for instance, £jv = perturbation theory, 



with 



5e 



(2) 



E 



Sk - £7 



(40) 
(41) 



When going back to the inverse compressibility, one 
should bear in mind, however, that the variation of VeS 
when going from iV— 1 to N electron is the same as when 
going from to A^-|-l. As a consequence, it can be 
checked easily that the linear variation Sef'^ cancels for 



At this level of approximation, the variation of the po- 
tential inside the dot is just the constant S^^, which will 
just shift the one particle energies but not give rise 
to any fluctuation. One should keep in mind though 
that Eq. ( |33|) applies only at places where nGTF(r) 7^ 0. 
Therefore the boundaries of the effective potential can 
be affected, in a way that depends, of course, largely on 
the external confining potential, and in particular on the 
second derivative of this latter at the the boundary. If 
the curvature is weak, fairly large displacements of the 
boundary can occur, which could play a significant role. 

If, on the other hand, one assumes a billiard-like con- 
fining potential, the only source of modification of the 
effective potential comes from (Stt-gtfi which is obtained 
through the equation 

J dr'^(r,r')<5n°GTF(r')+ / dr'^(r, rOfc^^^F = 



all levels except the two top ones, 
proximation one therefore gets 



At this level of ap- 



N 



(42) 



where A^E"^^ is the constant interaction model result 



given in Eq. (||), and 



N-l 



^scramb 



Sen+i +SeN+Y^ ^'^^3 



is the correction due to scrambling. The levels used for 
A^Ej^^ are those for the TV-electron Vcs- Throughout 
this scrambling argument we assume that the filling of 
the levels does not change. 



First order terms 



(37) 



Noting that J {SVcs / 6n)SnQj,p is, up to an inessential 
additive constant, the variation SVcS of the effective po- 
tential we are interested in, and that 



The mean of the first order terms is zero by construc- 
tion. The variance of (Se,-^^ is 



var 



= A^ /drdr'|(i^,(r)?/.*(r')>l''5^(r)'5n(r') 



(38) 



we get 



(44) 

Using the wavefunction correlations ( [T^ ) and approxi- 
mating J'^{kp\v — r'l) ~ l/{-nkp\r — r'|) plus oscillating 
terms (fcj^lr — r'| > 1), we find 



5V,s{T)^-A{6n{Y)-5n) 



(39) 



where 6h = A5n^r^-p(Y) is a smooth function of order 1 
[A is the area of the dot). Note that we choose SVcn with 
mean value zero. This ensures that there is no change in 
the mean energy levels; any such change in the mean 
should be incorporated in the charging energy. 



a ■ 



(45) 



The prefactor depends on the geometry of the dot under 
consideration: 



= / d^nd^u' 



1 



— 5h{u)5h{vL) 



(46) 
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is a dimensionless quantity of order 1 (note vl = t/L). For 
a circular dot, for instance, with diffuse boundary scat- 
tering to make it chaotic, 5n = {[1 — (r/i?)^]^^/^ — 2}/2 
where R is the radius; in this case, a ~ 1.2. 

Bearing in mind that kpL c::^ (taking into ac- 

count the spin degeneracy), we see that the typical am- 
plitude of the variation of the 5ei scales as N~^/'^ times 
the mean level spacing. This may, in fact, be numeri- 
cally significant for experimental dots — for a circular dot 
with 100 electrons, for instance, the magnitude would be 
~ 0.12A . 



Second order terms 

The variance of the sum 'Y^^ ^ can be evaluated 
along the same lines. Since it turns out to be parametri- 
cally smaller, we will be less careful with the constant. 

First we note that the fluctuation of an individual ma- 
trix element is 



m5VMj\)) 



dr c?r ' Jo ( fci (5r ) Jo ( fc j (5r ) (5n ( r ) (5n (r ' ] 

(47) 



where 5y = r— r'. Following the same reasoning as above, 
this implies that 







if SkL > 1 



aA'^/nkFL if SkL < 1 



(48) 



where 6k — ki — kj. 

This result is now inserted in the expression m\) for the 

(2) . — 

energy shift Se^ and then for the inverse compressibility 
Eq. (^3|). Assuming the matrix elements (j|^Vcff|j|) are 
Gaussian distributed independent variables, we find 



(49) 



Because of the saturation of the fluctuations of the levels, 
point (iii) in the random matrix criteria of Section the 
sum over levels is cut off at Nc — kpL. The final result is 




(Tn cx In kf.L 



A 

kpL 



(50) 




FIG. 3. The probability density of CB peak spacings including the scrambling effect on the top two levels as well as the 
fluctuations in the residual interactions. The total distribution (solid) as well as that for A*' even (dashed) and odd (dotted) 
is given for two sizes of quantum dot (A*' = 200, left column, or 80, right column) and two strengths of interaction {rs — l, top 
row, or the case of a zero-range interaction, bottom row). Note that the clear odd/even effect for Ts = 1 and the relatively small 
A'' = 80 remains visible even for the stronger zero-range interaction. The spacing is in units of the mean level separation A, and 
the origin corresponds to the classical spacing e^/C. 
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which is parametrically smaller than the first order re- 
sult (|5|). We therefore see that the the effect of scram- 
bling is dominated by the variation of the top levels, and 
that this adds contributions to the peak spacing fluctu- 
ations which while parametrically small on the scale of 
A could be numerically important for dots which are not 
too large. 



C. Numerics 

It is a simple matter to incorporate the scrambling of 
the top two levels into our numerical investigation of the 
ground state energy specified by Eq. (|^) . In random ma- 
trix theory, the change in an energy level upon varying 
a parameter is uncorrelated with its value [ |6T|j62| , so we 
simply choose random Ssn and Sen+i from a Gaussian 
distribution with variance Eq. (^5|). The ground state 
configuration is then found as before. 

Fig. 3 shows the resulting probability density of CB 
peak spacings. In the Vs = 1 case there is still a sub- 
stantial difference between the distributions of even and 
odd spacings: for N even, the large spacing tail com- 
ing from the level-spacing distribution (Wigner surmise) 
survives even for less than 100 electrons. When the inter- 
actions are stronger, as in the zero-range interaction case 
shown, the even and odd distributions are more similar, 
but there is still a discernible difference. The magni- 
tude of the fluctuations increases further: we now find 
rms{A^EN)/A = 0.33-0.32 for r, = l and TV = 80-200. 

VI. BEYOND THE GAUSSIAN MODEL 

Our treatment of the fluctuations of the residual inter- 
actions in the previous sections relies on a certain number 
of assumptions. We have for instance neglected correla- 
tions among the Afy and Ni'ji (except for the obvious 
Mij = Nij for a zero-range interaction) when we know, 
for instance, that the correlation coefficient of with 
Nij is of order 1 even in the infinite size limit. Similarly 
we have assumed the absence of correlation between these 
residual interaction terms and the one particle energies e^. 
Furthermore, we have used expressions such as Eqs. (p3|)- 
(^) which are derived in the large kpL limit. To check 
that these assumptions do not drastically modify the re- 
sulting distribution of CB peak spacings, we shall in this 
section implement the random matrix ensemble described 
in Section || [paragraph containing Eq. (|lj)] by directly 
performing numerical Monte Carlo calculations. 

To implement the constraints of our random matrix 
model, we start by considering the system to be a square 
of size L, for which we can use a plane wave basis 

Wx^ny) = exp[i27rri2;a;/L] exp[i27rrij,t//L] (51) 



ordered by increasing momentum. In this basis, we write 
the one particle Hamiltonian as H = D + V where D 
is diagonal with Da — i. V is taken to be a banded 
random Hamiltonian with independent complex matrix 
elements which are Gaussian distributed with variance 
vf^ — a{i)exp{~b{i)\i — For a{i) in the correct range 
p3[ such an ensemble yields (1) eigenvalues with a mean 
spacing close to one, (2) local fluctuations which are typ- 
ical of the GUE, and (3) eigenfunctions which are local- 
ized in energy. The localization in energy can be char- 
acterized by the inverse participation ratio in the plane 
wave basis: Ij = X^n n \{^x,ny\^j)\'^- Our requirement 
Eq. ( |l^ ) implies the constraint 

= m){k,L) . (52) 

In practice, we use a simple algorithm to tune the value 
of a{i) and b{i) in such a way that this constraint holds. 
Diagonalizing the random matrix specifled leads to eigen- 
values and eigenvectors which are then used to evaluate 
ground state energies for different number of particles via 
Eq. (^, and hence arrive at the CB peak spacings in this 
model. 

In this approach, deviations from the wavefunction cor- 
relations (|l^) required by normalization are naturally 
present. Furthermore, correlation between eigenvalues 
and eigenfunctions is present as well. Of course, in a 
real physical system, these effects would be automati- 
cally included because of the scattering at the boundary 
of the system — the transit time r, for instance, gives both 
the scale over which fluctuations saturate (Ti/t) and the 
length scale at which the Bessel function correlations be- 
come suspect (vft). We believe that our random matrix 
model includes the main features of these effects in a 
generic way. 

Fig. 4 shows numerical results for this model. We 
consider only the smaller value oi N, N = 80, because 
the effects are expected to be only significant for smaller 
quantum dots. On the other hand, we present results 
both with and without the scrambling effect of Section 
, and for different interaction strengths — Vs = 1 corre- 
sponds to a typical value for a GaAs quantum dot while 
the zero-range interaction case corresponds to the max- 
imum effect possible within a perturbative approach to 
the interactions. When these results are compared to 
those of Figs. 2 and 3, we see that there is a small 
but noticeable effect. This is highlighted in the quantile- 
quantile plots in the insets to Fig. 4 where an ordered 
set of RMT spacings is plotted against a similar set for 
the Gaussian model. The present RMT yields a more 
pronounced large spacing tail than the Gaussian model. 

VII. CONCLUSIONS 

In this paper we have studied the Coulomb blockade 
peak spacings for large but not infinite quantum dots 
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FIG. 4. The probability density of CB peak spacings within a random matrix model which goes beyond the Gaussian 
approximation. For A'^~80, the total distribution (solid) as well as that for A'^ even (dashed) and odd (dotted) is given for two 
strengths of interaction (rs = 1, top row, or the case of a zero-range interaction, bottom row) in the presence of scrambling (right 
column) and without scrambling (left column). Features of the distribution are more pronounced than in the Gaussian model 
(Figs. 2 and 3), suggesting that smoothing effects are slightly overestimated in the Gaussian approximation. Insets: Quan- 
tile-quantile plots of the random matrix model versus the Gaussian model; note the differences at large spacings. Throughout, 
the spacing is in units of the mean level separation A, and the origin corresponds to the classical spacing e^/C 



in the case that the interactions can be treated pertur- 
batively and the single-particle properties described by 
random matrix theory. Motivated by the Strutinsky ap- 
proach to finite Fermi systems ]50|-|5^ , we use an expres- 
sion for the ground state energy [Eq. (^] which includes 
the charging energy, the single particle quantization ef- 
fects, and the residual (screened) interactions treated in 
first order (Hartree-Fock). 

We presented results in four cases. First, in the limit 
kpL — !■ oo, fluctuations of the residual interactions are 
negligible. The CB peak spacing distribution has a dis- 
continuity at e^/C and a (5-function because of the odd 
N case (Fig. 1). The distribution is very different from 
that predicted by the simple CI model which neglects 
some terms of order A while keeping others. 

Second, fluctuations in the residual interactions con- 
tribute a term typically of order A/^/kpL to the CB peak 
spacings. The sharp features of the infinite dot limit are 
rounded out (Fig. 2). 

Third, the addition of an electron in the CB process 
changes the confinement potential in the quantum dot 
and so changes the single particle properties. This scram- 
bling effect also contributes a term typically of order 



A/\/kpL to the CB peak spacings, and the distribution is 
rounded further. However, the difference between N even 
and N odd remains substantial at large spacings where 
a tail is present for N even (Fig. 3) which is a remnant 
of the Wigner surmise distribution for level spacings. 

Finally, we took a first step towards investigating cor- 
rections to the Gaussian uncorrelated model for the fluc- 
tuations of the residual interactions. We found that the 
large spacing part of the distribution is somewhat under- 
estimated in the Gaussian model. 

Before ending this paper, we compare our theoretical 
distributions to those obtained experimentally [p5[-p9|. 
The observed distribution has little structure, being 
Gaussian near its peak but with longer tails for both 
larger and smaller spacings. The width of the distribu- 
tion is contentious: early work claimed that the width 
scaled with the charging energy j2^j2^ but this has not 
been verified by more recent experiments p7|-|2^ . In the 
more recent experiments, the width is clearly related to 
the mean level separation A and in the range 0.5-1.5A, 
somewhat larger than given by, e.g., the CI model. 

The theoretical results presented here certainly look 
more like the experiments than the benchmark CI model: 
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the distribution is more rounded and so less structured, 
and the odd/even alternation is much less pronounced. 
We suspect that further development of this approach 
will produce qualitative agreement with the experiments. 

In the end, however, the theory presented here does 
not provide an adequate quantitative description of the 
experiments. To start with, the magnitude of the fluctu- 
ations is actually less than in the CI model (0.58 versus 
0.32 for TV = 200) and so differs from experiment. With 
regard to the odd/even effect, the situation is somewhat 
less clear, and so more interesting. Looking at Fig. 3, for 
instance, one observes that a substantial odd/even effect 
is still present for = 1, while it has been largely washed 
out in the zero-range limit of the screened interaction. It 
should be stressed that this drastic modification of the 
distributions can be traced back mainly to the relatively 
modest change in the mean values (Mu) — A/2 and {Nij), 
Eqs. (Ill)-®, which change from about 0.3A for = 1 
to 0.5A in the zero-range limit. The sensitivity to the nu- 
merical value of these quantities should not come as a sur- 
prise since we know that the Stoner instability, for which 
a spontaneous polarization of the spin of many electrons 



takes place, occurs for {Ni, 



A. These observations 



suggest an experiment: placing a good conductor close 
to the 2DEG in such a way that it can screen the inter- 
action between electrons within the quantum dot should 
significantly affect the peak spacing distribution. If the 
auxiliary screening is strong enough, a strong odd/even 
effect should be visible. A double quantum well structure 
seems like a possible setting for such an experiment. 

We can divide the potential explanations for the quan- 
titative discrepancy between theory and experiment into 
basically three categories. The first one, which we can- 
not absolutely exclude, is that noise in the experiments 
is corrupting the results p7| , |64| . The second is that ef- 
fects not included here, but interpretable within the kind 
of Fermi liquid theory on which we have relied through- 
out this paper, dominate. Some candidates are (a) the 
effect of the changing gate potential on the single parti- 
cle levels [|l^,|5| — this would, however, somewhat contra- 
dict the remarkable stability observed in the wavefunc- 
tion — (b) boundary effects in the scrambling which 
need to be considered using a realistic description of the 
dot, and (c) consequences of non-chaotic dynamics within 
the dots |ll|,|l|]. Finally, a third completely different 
kind of interpretation is that the perturbative approach 
to interactions taken here is not accurate even though 
rs~ 1-2 p5| , pT| , ^j4Tl . Further experimental and theoret- 
ical work is needed to discern between these options. 
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APPENDIX A: VARIANCE OF THE Mij 

In this appendix we sketch the argument which leads 



to Eqs. (23)-(E6|). As in Section VI, we assume that the 



configuration space is a square of size L, for which the 
plane wave basis |k) — \nx,ny) introduced in Eq. ( [5l| ) 
can be used. 

Expanding the eigenfunctions in this basis. 



k 



(Al) 



we can write the My 's as 



^'^i] «iki<k2''^jk3Wjk4'^ki-k2,-k3+k4V'TF(ki - kz) 

ki,k2,k3,k4 

^J2^tf{cl)W,^W*^ (A2) 



where we define 



ki-k2=q 

Our random matrix model (Section ||) implies 



(A3) 



{vtkiVjkJ = {kpL) ^*y(5kik2 if 5k < n/L (A4) 
= ff (5fc > n/L 



with Sk = ||ki| — ki\ and ki ~ ^/2/mEi/h. From this, we 
deduce 



var(M.^,) =^ ^ E ^TF(q)(|W^^ql')(|W^.ql 



(A5) 



|Wjqp can be interpreted as (27rfci)~^ times the area 
of the intersection of two rings of diameter ki and width 
2'k/L, centered at a distance |q| — q. Simple geometry 
therefore gives, for 2tt/L < |q| < 2ki — 2tt/L 



{qL)^{2hLY - [qLY 
We obtain for i ~ j 



(A6) 



The variance of Mu and Nij and the covariance be- 
tween Mij and Nij can be computed along the same lines. 

[1] Single Charge Tunneling: Coulomb Blockade Phenomena 
in Nano structures, edited by H. Grabert and M. Devoret 
(Plenum, New York, 1992). 



14 



[2] Mesoscopic Electron Transport, edited by L. L. Sohn, 
L. P. Kouwcnhoven, and G. Schon (Kluwor, New York, 
1997). 

[3] D. K. Ferry and S. M. Goodnick, Transport in Nanostruc- 
tures (Cambridge University Press, New York, 1997). 

[4] Nanotechnology, edited by G. Timp 

(Springer- Verlag/AIP, New York, 1999). 

[5] M. A. Kastner, Rev. Mod. Phys. 64, 849 (1992). 

[6] L. P. Kouwenhoven, C. M. Marcus, P. L. McEuen, S. 
Tarucha, R. M. Westervelt, and N. S. Wingreen, in Meso- 
scopic Electron Transport, edited by L. L. Sohn, L. P. 
Kouwenhoven, and G. Schon (Kluwer, New York, 1997), 
pp. 105-214. 

[7] D. C. Ralph, C. T. Black, J. M. Hergenrother, J. G. 

Lu, and M. Tinkham, in Mesoscopic Electron Transport, 

edited by L. L. Sohn, L. P. Kouwenhoven, and G. Schon 

(Kluwer, New York, 1997), pp. 447-467. 
[8] I. L. Aleiner, P. W. Brouwer, and L. I. Glazman, in prepa- 
ration (unpubHshed). 
[9] R. A. Jalabort, A. D. Stone, and Y. Alhassid, Phys. Rev. 

Lett. 68, 3468 (1992). 
[10] M. Stopa, Phys. Rev. B 54, 13767 (1996). 
[11] M. Stopa, Physica B 251, 228 (1998). 
[12] G. Hackenbroich, W. D. Heiss, and H. A. Weidenmiiller, 

Phys. Rev. Lett. 79, 127 (1997). 
[13] Y. Alhassid, M. Gokgedag, and A. D. Stone, Phys. Rev. 

B 58, 7524 (1998). 
[14] R. O. Vallejos, C. H. Lewenkopf, and E. R. Mucciolo, 

Phys. Rev. Lett. 81, 677 (1998). 
[15] R. O. Vallejos, C. H. Lewenkopf, and E. R. Mucciolo, 

Phys. Rev. B 60, 13682 (1999). 
[16] E. E. Narimanov, N. R. Cerruti, H. U. Baranger, and S. 

Tomsovic, Phys. Rev. Lett. 83, 2640 (1999). 
[17] L. Kaplan, Phys. Rev. E 62, 3476 (2000), 

[arXiv:nUn.CD /0003013] . 
[18] E. E. Narimanov, H. U. Baranger, N. R. Cerruti, and 

S. Tomsovic, submitted to Phys. Rev. B, arXivxond- 

mat/0101034 (unpublished). 
[19] S. R. Patcl, D. R. Stewart, C. M. Marcus, M. Gokcedag, 

Y. Alhassid, A. D. Stone, C. I. Duruos, and J. J. S. Har- 
ris, Phys. Rev. Lett. 81, 5900 (1998). 
[20] A. M. Chang, H. U. Baranger, L. N. Pfciffcr, K. W. West, 

and T. Y. Chang, Phys. Rev. Lett. 76, 1695 (1996). 
[21] J. A. Folk, S. R. Patel, S. F. Godijn, A. G. Huibers, 

S. M. Cronenwett, and C. M. Marcus, Phys. Rev. Lett. 

76, 1699 (1996). 
[22] M. G. Mayer, Phys. Rev. 60, 184 (1941). 
[23] R. Latter, Phys. Rev. 99, 510 (1955). 
[24] O. Bohigas, X. Campi, H. Krivine, and J. Treiner, Phys. 

Lett. 64, 381 (1976). 
[25] U. Sivan, R. Berkovits, Y. Aloni, O. Prus, A. Auerbach, 

and G. Ben-Yoseph, Phys. Rev. Lett. 77, 1123 (1996). 
[26] F. Simmel, T. Heinzel, and D. A. Wharam, Europhys. 

Lett. 38, 123 (1997). 
[27] S. R. Patel, S. M. Cronenwett, D. R. Stewart, A. G. 

Huibers, C. M. Marcus, C. I. Duruz, J. S. Harris, K. 

Campman, and A. C. Gossard, Phys. Rev. Lett. 80, 4522 

(1998) . 

[28] F. Simmel, D. Abusch-Magder, D. A. Wharam, M. A. 
Kastner, and J. P. Kotthaus, Phys. Rev. B 59, R10441 

(1999) . 



[29] S. Liischer, T. Heinzel, K. Ensslin, W. Wegscheider, and 

M. Bichler, arXiv;cond-mat/0002226 (unpublished). 
[30] T. T. Ong, H. U. Baranger, C. M. Marcus, and S. R. 

Patel (unpublished). 
[31] R. Berkovits, Phys. Rev. Lett. 81, 2128 (1998). 
[32] P. Walker, G. Montambaux, and Y. Gefen, Phys. Rev. 

Lett. 82, 5329 (1999). 
[33] P. Walker, G. Montambaux, and Y. Gefen, Phys. Rev. B 

60, 2541 (1999). 
[34] A. Cohen, K. Richter, and R. Berkovits, Phys. Rev. B 

60, 2536 (1999). 
[35] K.-H. Ahn, K. Richter, and I.-H. Lee, Phys. Rev. Lett. 

83, 4144 (1999). 
[36] N. W. Ashcroft and N. D. Mermin, Solid State Physics 

(Holt, Rinehart and Winston, New York, 1976), p. 344. 
[37] M. Koskinen, M. Manriinen, and S. M. Reimann, Phys. 

Rev. Lett. 79, 1389 (1997). 
[38] I.-H. Lee, V. Rao, R. M. Martin, and J.-P. Leburton, 

Phys. Rev. B 57, 9035 (1998). 
[39] K. Hirosc and N. S. Wingreen, Phys. Rev. B 59, 4604 

(1999) . 

[40] O. Prus, A. Auerbach, Y. Aloni, U. Sivan, and R. 

Berkovits, Phys. Rev. B 54, R14289 (1996). 
[41] R. Berkovits and B. Altshuler, Phys. Rev. B 55, 5297 

(1997) . 

[42] Y. M. Blaritcr, A. D. Mirlin, and B. A. Muzykantskii, 

Phys. Rev. Lett. 78, 2449 (1997). 
[43] P. W. Brouwer, Y. Oreg, and B. I. Halperin, Phys. Rev. 

B 60, R13977 (1999). 
[44] H. U. Baranger, D. Ullmo, and L. 1. Glazman, Phys. Rev. 

B 61, R2425 (2000). 
[45] I. L. Kurland, I. L. Aleiner, and B. L. Altshuler, Phys. 

Rev. B 62, 14886 (2000). 
[46] P. Jacquod and A. D. Stone, Phys. Rev. Lett. 84, 3938 

(2000) . 

[47] P. Jacquod and A. D. Stone, cond-mat/0102029 (unpub- 
lished) . 

[48] M. Gutzwiller, Chaos in Classical and Quantum Mechan- 
ics (Spring- Verlag, New York, 1991). 

[49] B. L. Altshuler and B. D. Simons, in Mesoscopic Quan- 
tum Physics, edited by E. Akkermans, G. Montambaux, 
J.-L. Pichard, and J. Zinn-Justin (Elsevier Science, Am- 
sterdam, 1995), pp. 1-98. 

[50] V. M. Strutiriskii, Nucl. Phys A 122, 1 (1968). 

[51] M. Brack, J. Damg0ard, A. S. Jensen, H. C. Pauli, and 
V. M. Strutinsky, Rev. Mod. Phys. 44, 320 (1972). 

[52] D. Ullmo, T. Nagano, S. Tomsovic, and H. U. Baranger, 
arXiv:cond-mat /0007330 (unpublished) . 

[53] Chaos and Quantum Physics, edited by M.-J. Giannoni, 
A. Voros, and J. Zinn-Justin (North-Holland, New York, 
1991). 

[54] M. V. Berry, J. Phys. A 10, 2083 (1977). 

[55] M. Srednicki, Phys. Rev. E 54, 954 (1996). 

[56] O. Agam, N. S. Wingreen, B. L. Altshuler, D. C. Ralph, 

and M. Tinkham, Phys. Rev. Lett. 78, 1956 (1997). 
[57] I. L. Aleiner and L. I. Glazman, Phys. Rev. B 57, 9608 

(1998) . 

[58] Y. M. Blanter, Phys. Rev. B 54, 12807 (1996). 
[59] Y. M. Blanter and A. D. Mirlin, Phys. Rev. B 55, 6514 
(1997). 

[60] P. Leboeuf and A. Monastra, Phys. Rev. B 62, 12617 



15 



(2000). 

[61] M. Wilkinson, J. Phys. A 22, 2795 (1989). 
[62] E. J. Austin and M. Wilkinson, Nonlinearity 5, 1137 
(1992). 

[63] Y. V. Fyodorov and A. D. Mirlin, Phys. Rev. Lett. 67, 



2405 (1996). 

[64] T. Heinzel, private communication (unpublished). 
[65] D. R. Stewart, D. Sprinzak, C. M. Marcus, C. I. Duruoz, 
and J. S. Harris, Science 278, 1784 (1997). 



16 



